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ABSTRACT 


The small signal dynamic stability of power systems 
has been a subject of interest for some twenty years. 
It continues to grow in importance as the control 
requirements of the generating stations become more 
sophisticated and demanding. In this thesis, the 
small signal performance of a regulated synchronous 
machine connected to an infinite bus is described by 
a set of differential equations of the form (x) = 

(A) (x). The transformation of the system into 
Schwarz form is used and the st:ability conditions for 
the system under study are established. The closed 
stability regions are obtained as the result of this 
study and the effects of various system parameters on 


stability regions are studied. 
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INTRODUCTION 


1.1 Background 


Electrical power systems, considered from the point of view of their 
electro-mechanical operation, are very complicated. Power system 
engineers have devoted much thought and effort to stability studies 
since about 1925. The tendency of a power system or its component 
parts to develop forces to maintain synchronism and equilibrium is 
known as stability. In general, stability studies are classified by 
whether they involve steady-state or transient conditions. Dynamic 
stability is the term associated with the small signal performance 
and is applied to operations above the ordinary steady-state limits. 
It can be realized by the use of automatic control devices such as 
voltage and speed regulators. The small signal dynamic stability of 
electric power systems has been a subject of major theoretical and 
practical interest for some twenty years. It continues to grow in 
importance as the control requirements of the power plants become 
more sophisticated and demanding. 

Power system stability studies as usually conducted, with all of their 
assumptions and approximations, are far from an exact science. 
Possibly some of the approximations can be justified because of the 
power system conditions analyzed. However, they will almost never 
correspond to the actual conditions existing on the system when the 
actual disturbance occurs. Also, it is virtually impossible to know 


how the actual power system loads vary with frequency and voltage 
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magnitude. It is certain that all the loads do not act as fixed 
impedances. It is equally certain that they do not remain as constant 
current loads. Possibly, sufficient information will never be avail- 
able on loads to represent them, and even if the information were 
available, it might be too complex to represent them accurately in a 
practical digital study. Usually it will not be possible to represent 
all of the generating plants connected to an interconnected system. 
This results in further approximation since there seems to be no exact 
simple equivalent that represents the generating plants. Further, the 
representation of individual generating plants also introduces errors 
because the machine impedances are not really constants. It appears 
that this can result in a paradoxical situation where the system will 
be calculated to be transiently stable and also calculated to be un- 
stable under steady-state. Most analytical treatments [9, 10] of 
stability problems have been restricted to one or two generating sets 
because of the computational difficulties involved in applying the 
older methods. However, the techniques of modern control theory have 
partly removed this difficulty, subject to the requirement that the 
system be described by a set of differential equations in the state- 
space form 

[x] = [A] [x] (1) 
Laughton [8] has proposed a method of obtaining the [A] matrix of a 
multi-machine power system by using matrix elimination techniques to 
extract [A] from the complete algebraic and differential equations of 
the whole system. The paper by Undrill [7] describes a method of 


building up the [A] matrix of the multi~machine power system from the 


3 ata : | 
| Bead 26, 405 Jon ob abso Site hp ats ntetieo 2b aT -ebusitngam 
Jnatenes as stent Jon ob ysl) dad mbgatae vitaupe vat 3T wee ; 


etsw sokismrotnt efj 2! save bre ,wod3 Imsaevges 03 abgol go elds 

b wt ylerstuoss mani 10sa07q81 64 xslqmas dos sd trighm 32 ~eidsliava 
jnagetge: of sidtezog 5d.3jon [Ibw 3t yilevall -xbuse Letigth Isotiossq 
«eiaye beJosnaozrs3ak ns 03 bsioevanoa ednslq antzerensg ala 20 Lis 
Jnex9 on od ot atlases sisiy sonte noltsmixorgqas. Tojxu? al etlveaex eldT 


x 
od3 ,xeditut |.esosiq anitetons, sii sinsasrqs+ defi Jnslaviupes olgmte . 
eter7s aon ita oals etnoaiq gntontonag isubivibri jo colisinsesrges 
strneqqre tI ‘qaiabireno uifsst jon ots egonbbagut sntdsam ait saunoed 
fitw mateye sd4 siadw Aoksauats fsoixobsisq 5 oi Yiuesi oso abe ter3 a 
“fi 90 02 betalustag cele bas sidese ylanotenst2 sd o2 bsdeduateo ad : 
io [LOL .2} adnamgssrt [sottylans secM -23gj2-ybsete tabmy sidase 
2398 gidtiatensg dvi 10 s00 03 fesotytesy nasd syed emoddorg Yatiidese : 
stiz ghivtaae ak beviovnt eaksivoti?ls Isnottayiqmon ont to sauncsd 7 
avad yroid3 fotinos mrsbom to esupiongss sit ,ysvewoH .abodzer zebfo rs 
od3 jsd3 insmeziupst sda o3 Jostdbe .ytlwotttib atds bevomsy yistsq - _ 


"93838 Sj ni anoltsups fakinsrsi?ib to jae s yd Sediaoasb od asveye 
@201 s7aq8 

@ [x} Lal = (3) 
& to #itgen [A] sda gotatsade to bodjem.s besegetq ead [8] notdguad 
09 asupintioss aobisnimzia xiijem gate yd ootae x9wog setdsem-hitems 
Yo aohieupe Istsnes9I2ib ‘bas 2axduphy esplqoo so moyd. [A] toeazae 
20 hoddem & aedingads (S] ILiabay yd soqgeq Su” | adage sodw ade 
MP mak ete sone ‘entitonn- tn aso akedha. 14) aaa 
De eee i " ty oe ' . a 


-) 
‘ 


? 4 
aa | 


submatrices describing individual elements of the synchronous system. 
The matrix building approach offers significant savings in computer 
storage in comparison with the matrix elimination approach of 
Laughton [8]. The transient response of a dynamical system for small 
perturbations about an equilibrium (operating point) is completely 
described by a set of such equations. In addition, using the digital 
computer it is possible to decrease the time necessary for computation 
of stability analyses in this form. 

In this thesis, the dynamic system which consists of a regulated 
synchronous machine connected to an infinite bus is taken for the 
purpose of the study. The system is expressed in the vector-matrix 
form (1) and the characteristic polynomial is obtained. The trans- 
formation into Schwarz form is performed and the necessary and suffi- 
cient conditions for stability established. Using these conditions, 
the stability regions are obtained and the change of these regions 


with variation of system parameters is observed. 


1.2 Objective of Thesis 


The stability study of the individual unit can be carried out and the 
necessary and sufficient conditions for stability established, 
assuming the rest of the system as an infinite bus. These can be 
used to find the stable regions for the system under study. 
Liapunov's criteria for stability and eigenvalue methods presented in 
the literature [13, 14] give satisfactory results for smaller systems 
but become too cumbersome for large systems. It seems reasonable 


therefore to look for conditions which could simplify the matrix 
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representation of a system, i.e. minimize the computational effort. 


Consequently, the following objectives can be formulated. 


1. To find the necessary and sufficient conditions for stability 
2. To find the stability regions for the given system 
3. To look for the possibility of extending the procedures employed 


in 1. and 2. to multi-machine systems. 


Generally, the necessary and sufficient conditions for stability 
should yield the closed stability regions which should be of parti- 
cular interest since no paper has been published so far on this sub- 


ject. 
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CHAPTER iL 


MATHEMATICAL MODEL 


ote. Description of the System 


The dynamic system taken for the purpose of this study consists of a 
regulated synchronous machine connected to an infinite bus. The 
system is generally considered to be non-linear and can be expressed 


in the vector-matrix form as follows: 
x = £(x) (2) 


The individual generating unit connected to an infinite bus has been 
investigated by Yu and Vongsuriya in [3] ana by Kasturi and Doraraju 
in [4] in the configuration as shown in Fig. 1. The system under 
study consists of the synchronous machine, the voltage regulator, the 
speed govermor and the tie-line. In addition, the machine with its 
control equipment is connected to an infinite bus. The steady-state 


phasor diagram for the system under study is shown in Fig. 2, where 


eg - armature voltage in d axis 

re - armature voltage in q axis 

Qe, - armature terminal voltage 

e - infinite bus-bar voltage 

ey field winding applied voltage 

6 - load angle 

r - tie-line resistance between generator and infinite bus 


xX - tie-line reactance between generator and infinite bus 
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Figure 2 
Xq - synchronous reactance in d axis 
Xa - synchronous reactance in q axis 
I - the phase current 


The mathematical equations describing the state of the model at any 


instant consists of 


A. The control system equations 
B. Power transfer equations relating the mechanical input and 
electrical output power 


C. Machine equations 


In writing the equations governing the system under consideration, 


the following assumptions are considered: 
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1. Each stator winding is distributed to produce a sinusoidal 
magnetomotive force wave along the air gap 

2. Stator slots produce negligible variations in the rotor 
inductances 

3. The armature resistance of the synchronous machine is neglected 

4, Damping due to damper bars has been neglected 

5. Transformer voltages in Park's equations have been neglected 
compared to speed voltages 

6. Saturation in the machine is neglected 


7. The regulator is assumed to have no dead zone or limits 


The above assumptions are usual in power system studies and they are 


the same as those in [3] and [4]. 


242 ®e°The™ Control System Equations 


Assuming a conventional voltage regulator and speed governor, the 
control system can be written in a form which is commonly used in 


similar analyses [3, 4]: 
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where 
A - small change around initial operating point 
s - Laplace transform variable 


6 - instantaneous angular position of rotor 
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s6 —- angular velocity 
T, - mechanical power input to rotor 

= ~ 
us (Xoeq/T pe) He overall regulator gain 
T. - stabilizer time constant 
1 - exciter time constant 

— te 3 e ny 2. 

us ni Uedey overall stabilizer gain 
u - governor gain 


Tips oe governor time constants 


Uoyx 7 exciter gain 

Lol is stabilizer gain 

He = Masteety - regulator gain 

uy - convertor gain 

te 7 field winding resistance 

X ata’® akd* akq - mutual reactances between stator and rotor 


2.3 B. Tie-Line Voltage Equations 


The tie-line voltage equations derived from the steady-state phasor 


diagram shown in Fig. 2 of the system shown in Fig. 1 are as follows: 


eg = esin6d +r ig =X 
e = ecosdé+ri + Xi 
q q 
Ad pie GRP) p) 
e. eq oF ea 
Ou? 22 
i i ag a (5) 


Considering a constant frequency for the infinite bus, the instantan- 


eous angular displacement, velocity and acceleration of the machine 
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respectively are 


6 = 6 + 6 
s§ = s6 | + s6 = Wy + sé (velocity) 
s76@ = sty (acceleration) 


In these equations sO. ai is the synchronous speed and w is the 
undamped natural frequency. The system of equations (5) and the 
following system of machine equations can be obtained from [11], where 


also more detailed information can be found. 


2.4 C. Machine Equations 


Park's synchronous machine equations [11, 12] are used for the purpose 


of this study. The synchronous machine equations are as follows: 
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armature flux linkages in d and q axes, 
and damper flux linkages in d and q axes, 
respectively 
rotor self reactances 
armature and damper currents in d and q 


axes 


direct axis transient open-circuit time 


constant 


direct axis transient short-circuit time 


constant 


introduced here to complete the nomenclature 


of symbols used in the following text: 


ve - armature winding resistance in d 
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ay - electrical torque 
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or q axis 


P - real power output of machine 
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Q - reactive power output of machine 
H - inertia constant 

M - moment of inertia 

D - damping coefficient of machine 

W - undamped natural frequency 


In this study the stability of the power system due to small distur- 
bancies will be examined. The following nine equations are the usual 
equations for the similar stability studies obtained by linearizing 


the operational equations (3), (4), (5) and (6) about an operating 


point. 
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The quantities C40? eae Go. 6, etc. represent the steady-state 


to 


values at the particular operating point. Eliminating T,, are and 


e. in the preceding equations, five homogeneous equations are 


obtained, which can be put into the matrix form as follows: 


if 0 - Yr X -v cosé A Vq 
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h(s) = g(s) G(s) 
We 
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A characteristic equation of the form 
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can be obtained from the characteristic determinant of (8). Linear- 
ization of the system equations and the derivation of the character- 


istic equation is shown in more detail in [3] and [4]. 


vw tt We ae 
ne ‘aria 20 sonar es tno nto 79 a 194 
a, SMT ae fe? nx ays reps even sand 


= 
: : ’ 7 
te ed - 
_ o 
ter OF 
_ 


ai ei ea ue 
ai i. 
_ 7 7 : " 


i ES) 


CHAPTER 101 


CONDITIONS FOR STABILITY 


3.1L” Derivation of necessary and sufficient Conditions for Stability 


The theorems given by Ogata [1] will be stated. This theorem gives 
the answers to the problem of deriving necessary and sufficient 
conditions for stability of a linear dynamic system and will be em- 
ployed later to derive these conditions for the system under investi- 


gation. 


Theorem 8 - 9, [1] p. 466: 


Consider the linear time-invariant system 


y = Wy (11) 
where 
0 ] 0 0 ° 0 0 0 | 
| 
-b 0 1 (Os ices eas 0 0 0 | 
{ 
: | 
x; e 0 @) 
| 0 bo-1 0 ad 0 | 
| 
W = | ; (12) 
ee : 5 . 5 
| e e e | 
| 0 0) 0 0 aiaters -b, 0 1 
0 0 0 0 . 0 ah praia 
and bo > bos b, Py or ras bo are real quantities. The origin of the 


system is asymptotically stable if and only if 
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b, 0, b, E> Li b, >.0 seaistetete corer De ae a 


(W is called the Schwarz matrix) 


The above theorem gives directly the necessary and sufficient con- 
ditions for stability for the linear time-invariant system (11). The 
dynamic system under study is described by (1) and is of the form 

x = Ax. If we can find the similarity transformation and transform 
the system under study into the form given by (11), the whole problem 
of stability can be solved and necessary and sufficient conditions 

for stability established. In [2] we can find the confirmation of 
this assumption, since on the page 100 is stated as follows: 

"once the matrix A has been transformed into Schwarz form, the 
stability problem is solved immediately; the necessary and sufficient 
conditions for asymptotic stability being that all bis are positive 
(isle 3) nee on) es Another theorem from (1) will be dntroduced 
and the similarity transformation given by this theorem will be em- 


ployed later. 


Theorem 8 - 11, [1] p. 467: 


Ie @ Teeedl Ce@meceine melee (C 


C 1 Gis wheleern ene 0 | 

0 0 a eet «ohare ’ e's 6) | 

| , | 
CGa= | x a ° 
0 0 chatentieae rete 1 

a oH a | Sie elenel 6 one ay 4 
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is similar to the Schwarz matrix W. 
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c=T lyr (13) 


(T:non-singular), then the number of eigenvalues of C which have 
negative real parts is equal to the number of positive terms in the 


sequence 


pee babe gab aera oe one, bab 


Be Cle2 com a2. se 
provided that none of the bis is zero. 


Using this theorem, the following relationship will be used in order 


to transform the C matrix into the W matrix: 


w= tcT ! (14) 
where 
W - is Schwarz matrix 
C - Companion matrix associated with characteristic 
polynomial 
T - Transformation matrix 
Ral - Inverse of the Transformation matrix 


Generally speaking, every non-derogatory matrix A can be converted 
to the companion form and to the Schwarz form if the particular 


transformation matrix and its inverse can be found. 


The characteristic equation of the system under study (10) can be 


rewritten in the following form: 
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To outline further investigation at this point the following must be 
done in order to establish the necessary and sufficient conditions 


for the dynamic system: 


1. Find the companion matrix C of the system 

2. Find the transformation matrix T and its inverse 

3. Perform the matrix multiplication and obtain the Schwarz form 
4, Establish the necessary and sufficient conditions for stability 


for the system under study 


Using [1] and [2],the companion matrix associated with the character- 


istic polynomial given by (15) was found to be as follows: 


wey, 1 0 0 0 0 - 
, 0 1 0 0 0 o | 
| ° 0 0 1 0 0 " 
c= | 0 0 0 0 1 0 0 
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aC 0 0 0 0 1 t 
0 0 0 0 0 0 - 
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To obtain the transformation matrix T, two different methods can be 
followed. One is due to S. G. Loo and is shown in detail in [2]. 
The other method which seems to be more practical for the purpose of 
this study is that given by C. F. Chen and H. Chu and shown in [5]. 
These two authors also give the method for constructing the inverse 
of the transformation matrix [6] and this will be employed later in 


the text. 


Using the material [5], the transformation matrix T was found to be 


as follows: 


] 
1 0 0 0 0 0 c 
| 0 1 0 0 0 0 0 | 
{ 
rk | 
= 0 1 0 0 0 0 | 
Ns | 
| 
: | 
C | 
Sy i 
Tas 18 10 ae 0 A 0 0 0 | 
| 51 | 
| Cy,» 
| bat 0 = 0 1 0 0 
| Ch bd 
| 
c C 
| 0 me 0 see 0 1 0 
| C31 C31 | 
| : : | 
24 0 ae 0 = 0 1 | 
Coy Coy Coy | 


(17) 
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where che Fete = 1, 2, 3.......n/ are the elements of the Routh array 
for the given characteristic polynomial. The positions of the Routh 


array elements are indicated as matrix double subscript notations. 


As mentioned in the preceding text, 6 gives the material for ob- 
taining the inverse of the transformation matrix T; and it was found 


to be of the following form: 
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eee AEH mead 


2h 


where oie (i,j = 1, 2, 3.......n) are the same elements of the 
Routh array for a given characteristic polynomial as before. Further 
we have to construct the Routh array in order to evaluate the trans- 


formation matrix and its inverse. In [1] we can find the method for 


constructing such an array and it was found to be as follows: 


a! u a,/a, a,,/a, a,/a, 
r§ a,/a, a,/a, ac/a, ao/a, 
> ag b, b, 0 
at Cc, Cy C4 0 
ve d, d, 0 0 
ie a) e, 0 0 
r f, 0 0 0 
ae g, 0 0 0 
The coefficients for our seven-degree polynomial Ajo ays ag cere Ay 


are known and we have to express all the elements of this array in 

terms of these coefficients. In other words, we want to know b,, b,, 
u > — 

ba» Cis Cos Cas dj, dj, €)5 &5» £, and g, in terms of a,s Ge Oe 

2, 3 «eee» n). The following relationships hold. 
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by) 2 aa, aa3/aay 

ey ie “oii ii a,az/aia, 

b, “33 "ana. = a,a7/a ac 

ae Bs (aa, ~ aa) a, SAE lies a, ac) a,*/(a,a, - a.a,) aja, 
Co = (aa, - ajag) as* - (azag - aya) a,*/(aza, - a,a5) alas 
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OF ere Gy.coe = c,d) Gy 
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fy Serco. di e5/e, 

hate ee (19) 


These relationships allow the expression of all the elements of the 
Routh array in terms of coefficients of the characteristic polynomial. 
Now substituting the following values from the Routh array into T and 


t 1 matrices in order to get them to the suitable form for matrix 


multiplication (14): 


gh SIA% 8 C88 - gage) - “28 Cage ~ yee) 
oB\y8 


pa\go;d — Alps 


@\pa gd - edg> 


fB\ chro ~ 99 7b 
ofl 


. r2\o87b - gbys 


b. Gaia ps Seb Cea oie - 2 


Coy/Cyy = ar/ay C32) Cay = Pa! 
C33/C3, = b3/by Cya/Cyy = Co/cy 
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Ceo/Cg, = en/e, (20) 


Substituting equations given by (20) into (17), the matrix T is 
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Similarly, substitute the elements of the Routh array (20) into (18). 


To make the notation simpler, let 
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1 0 0 0 0 0 0 
0 1 0 0 0 0 0 
) 
alee 0 1 0 0 0 0 
M 
3 d, 
ud oes 0 cite 0 il 0 0 0 
1 
Co 
K 0 a 0 7 0 0 
Cc 
1 
by 
0 L 0 ——— 0 1 0 
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a 
-N 0 M 0 -= 0 1 
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Now matrices T and T ! are a suitable form for matrix multiplication 
and obtaining the Schwarz form (14) gives the relationship for trans- 
formation into Schwarz form and in this equation T is given by (21), 


T | by (22) and C by (16). After multiplication, the Schwarz matrix 


is 
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eect 0 1 0 0 0 0 
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When the system is transformed into the form given by (11), the 
necessary and sufficient conditions for asymptotic stability can be 
established. Theorem 8-9 introduced earlier in the text and given 

in [1] can now be applied in order to obtain the stability conditions 
for the dynamic system under study. It is stated in the above 
mentioned theorem that the necessary and sufficient conditions for 
asymptotic stability are that all bis are positive (i= 1, 2, 3 
..--e. n). If this is applied to the Schwarz matrix given by (23), 


the following seven conditions for stability can be found. 


12 e_/e rad) (24) 


The set of inequalities given by (24) forms the necessary and 
sufficient conditions for asymptotic stability for the dynamic 

system given by (15). It was mentioned previously that all constants 
involved in these conditions can be expressed in terms of coefficients 


of the characteristic polynomial, namely in terms of Ao» Ayo ges +87- 
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For the time being they may be left in the form given by (24) because 
this seems to be the simplest form in which to express the stability 
conditions for our case. Speaking generally, the stability conditions 
given by (24) must be satisfied simultaneously for some operating 
point and only then is it possible to assume that the dynamic system 
is stable at that particular operating point. Choosing some systems 
of coordinates, it will be possible to find stability regions in that 
plane if it is checked for stability at every point in that particular 
plane. In other words, for some point in the chosen plane to be 


stable the conditions given by (24) must be satisfied simultaneously. 


3.2 Multimachine Systems 


In the preceding text a single-machine system connected to aninfinite 
bus has been taken for analysis. However, with the state variable 
approach given by [7] and [8], the study can be extended to multi- 
machine systems. Speaking generally, the A matrix must be obtained 
for the considered system and the characteristic equation of the 
system established. Then the companion matrix of the investigated 
system must be found and using the transformation (14) the Schwarz 
form obtained. Theorem 8-9 can be applied again and the necessary 
and sufficient conditions for asymptotic stability of a dynamic system 
established. 

The method of analysis used in this study is quite general. To obtain 
the stability regions, the same technique can be used in multimachine 


systems as was used in the single-machine case. 
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Dynamic stability analysis of a large interconnected power system is 
extremely time-consuming and laborious and may even exceed the 
storage capacity of modern fast computers because of high order of 
the A matrix. Hence methods have been developed to obtain simplified 
models of the systems based on the speed of response of the variables 
or the nodes. One of such methods is given by Kappurajulu and 
Elangovan [15] and gives the reduced model of dynamic systems. It 
seems reasonable, therefore, to consider the possibility of reducing 
the system with more than three machines to some simplified form in 


order to simplify the stability calculations. 


3.3 Outline of the Method by Yu_ and Vongsuriya 


The steady-state stability limits of a regulated synchronous machine 
connected to an infinite bus were investigated by Y. N. Yu and 

K. Vongsuriya in [3]. The open stability regions were obtained as a 
result of this study and the method of the approach used in [3] will 
be outlined briefly in the following text. The system under investi- 
gation in this paper is that shown in Fig. 1 and the stability of the 
system due to small load disturbance is studied. Generally, the 
Routh-Hurwitz criterion and D-partition method are used in this 
steady-state stability study. The characteristic equation of the 
investigated system (10) is obtained and it can be written in such a 
form that the voltage regulator gain us and the stabilizer gain 


uy can be separated from the remainder of the equation. 
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To find the stability boundary on a Wi ue plane, let 


je ai ey 


which corresponds to the imaginary axis on the complex frequency 


plane. Separating P, Q and R in (a) into real and imaginary parts, 


P(ju) = P,(w) + 4P,(u) 


i] 


Q(jw) Q,(w) + 4Q,(w) 


R(jw) = R,(w) + 4R,(w) (c) 


the following can be obtained: 


me ee ue O70) Ra (a) = 9 0 
ut Pr 3) i ul Q,(w) + Rw) = 0 (d) 
Hence 
1 > R, (w) Q, (w) 
De ilar 
= R, Cw) Q, (w) 
; P, Cw) SOR CO)s 
2 Ey ee (e) 
where 
he ae ce 
P (w) Q, (w) 
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For a non-trivial solution to exist, 4 must be non-zero. Theoreti- 
cally, the complete stability boundary on the pee plane can be 
determined by varying w from zero to infinity. Practically, only 

a variation of w from zero to 0.02 is necessary to obtain all the 
useful information. A point test on stability by the Routh-Hurwitz 
criterion is necessary to determine which side of the boundary is 
stable. 

The effects of the saliency and the short-circuit ratio of the 
synchronous machine and that of the tie-line resistance and reactance 


on the stability of the system were investigated in this paper. 
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CHAPTER IV 


MODELLING OF THE SYSTEM 


eh System Parameters 


For the purpose of the study, the system will be introduced here with 
all the necessary machine data, regulator data and governor data. 
The system and its constants are the same as those used in [4]. The 


following constants are expressed in p.u. values and time in seconds. 


Machine data: 


P =1.0 M = .0337 

with 

PAS) Tho = 6.0 TA = lei3 
Kaas 2.0 Xq = 1.7 i. =) 6.0223 
x = 0.18 ene iat) 


Regulator settings: 


T = 1.12 T = 2.5 
e 


Governor settings: 


el eee iw = 10.0 


As stated earlier, a graphical display of the stable regions of the 
system under investigation is desired. For this purpose a Cartesian 
system, such that its X and Y coordinates correspond to us and ut 


which are the overall stabilizer and overall regulator gain 
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respectively, is used. These two parameters are controllable and 
were used in the similar study given by [3] to display the stability 


limits of a regulated synchronous machine. 


4.2 Initial Conditions of Synchronous Machine 


The initial values » v. and 6 
O ) 


*do? aot Vdo? Yqo° Yao? Yao? Vedo 
of a salient pole machine are found by the following expressions. 


In the steady-state, Park's equations become 


Vv = xi 


do q qo 
Vago = xad Vv XL 
Rp fdo d do 

E - waoedo Yqo qo 

Q a Ygo" do Vdo'ao 

See had v sino + tao x ta 
Ce = v,coss +e TY dao +L Ae 
ae = ve a8 eee 


and Veg are 
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P and Q are given. During the analyses in this 
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unknown, while veces 


thesis, P and Q will change but v, will be kept constant Chee = 


to 


1.0 = constant). The unknowns will be calculated from the following 
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relationships derived from the Park's equations for the steady-state. 


Pv 
‘to’ 


qo ea Comm ce go kee OF 
(P Shs + vie + Xo aye 


do q qo 
= et 2: 
iT “to Vdo 
a 2 
eeteepbcinthokee 
do Vv 
qo 
Ci apes ae i eee 
fd 
fo) qo d “do ee 
Yao = RY do 
Yao 3 Yqo 


6) = arctg 
Yes = x Loo =a tiG 


Further and more detailed information on deriving the initial con- 
ditions of a synchronous machine can be found in [3] and [4]. The 
initial conditions of the synchronous machine are calculated for all 


different power factors used in this study. These values are used 
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later in the process of establishing the coefficients of the 
characteristic polynomial, that is, they are used for the computation 


of the following constants: A A,» Sen 5 Hae Ags Ag: Results are 


il? 


tabulated in Appendix I. 


4.3 Coefficients for the Characteristic Polynomial 

The stability conditions were established in terms of the coefficients 
of the characteristic polynomial (24) of the dynamic system. The 

set of equations for calculating these coefficients is given here. 
From the characteristic equation (10) the coefficients of the poly- 


nomial are 
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gs a Th B. + i B, + i {B, + B, (T,, + rs ut) ys A, ee - 
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In the above set of equations, constants By, Bieeeeeeee Be, 


T at To Ty» To Ty and Ay, Ay seeeeeee Ag, Ay are introduced. 


Bie ) 


The equations specifying the above constants and the set of equations 
(28) are in agreement with [3] and [4]. 
To calculate the constants Bi» By ceseeees B. the following 


relationships can be used: 


) xl Ge oy ten) na 
B, = A, Xq + A, 
B, = A, Xq + Ay 
ooG 
Been eae <7 Tut Aven betes ya (Dao i==)ot Baad 
Bom (A. + AT!) y, + A, HA, TT. + By @ - — 
5 6 7a Xa 8 9 “do 2 W 


Ah 
eo 
Beige, Ay (D rat + A, (29) 


heh aes Be ‘ Ree ay 
eon ew wede te 


. es 

ks ilies - 

ie enna a ok 
sd nt ae hag ae ‘tne te ~ae 
0) 1 tae: 8 . 


, 


7 
i 
_ 

ce esa a 

a 


ee . ee | 
at yi a ee 
ee i | ats I ha a i 


= 
' 


dilehinlaneiens me 


- 


36 


where 


and 


for calculation of To Ty» T T, and T_, the following set 


od ae E F 


of relationships can be used: 


The 


The 


a L258 e 
Th = qT, + qT, T, au i ut 
ie = qT, T, + T, T; “ (T = T,) cre “ T, ut) 
Eee Saale T(t, fale ur) + (T, + Ly) Tale 
TE = qT, qT, + T, Ts + qT, Le 
Tp ee dads Ply 


last set of equations is for calculation of Ay» A, ceeee Ag» Ags 


following equations are used to calculate these constants. 
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The equations for the computation of initial conditions for syn- 
chronous machine were introduced earlier (27) and similarly the 
equations for the computation of all constants (31, 30, 29) necessary 


to evaluate the coefficients of the characteristic 
polynomial a), a5, a3 «e+. a 
As was mentioned before, it is convenient to leave these coefficients 


in terms of us and us (overall stabilizer and regulator gain 


respectively) and then to draw the stability regions for the 
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particular study on this plane. A digital computer program 
(Appendix II) is used to check as many points as required to 
establish stability regions. This program checks individual points 
on a grid. 

In summary, the following steps must be taken in order to obtain the 


stability regions for the dynamic system under investigation. 


1. Calculate the initial conditions of synchronous machine 


De 
Calculate constants To Th» To Ty Ty, and Ty 


3. Calculate constants Ay, Ap ceeeeeeee Ags Ay 


4, Calculate constants Bi» Bo, Ba, By, By and B, 


5. Establish the coefficients for the characteristic polynomial 


of the system, namely Aye Ay seeeeeeee ay 


6. Check for stability on WEue plane for as many points as is 


required. 
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CHAPTER AV 


RESULTS 


Stability Regions 


In this study the effects of the following parameters on the stability 


region of the system under study will be examined. 


6: 


Changes in power factor 

Changes in damping coefficient 
Stabilizer time constant effect 
Exciter time constant effect 

Effect of a one-time constant governor 


Effect of a two-time constant governer 


These studies are introduced in the following text and all the con- 


clusions drawn later in the text are based on these studies. 
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5.2 Effect of Different Power Factor on Stability Region 


In this study, the real power output of the machine is kept constant 
(P = 1.0 p.u.) in all investigated cases. The reactive power output 
of the machine is varied and the following distinct values of Q are 


used. 


The system data used here are those introduced earlier and given by 
(25). All constants are in p.u. values and time is in seconds. 
Using the set of equations (27), the initial conditions of the syn- 
chronous machine were calculated for all different values of Q and 
are shown in Appendix I in tabulated form. 

For each different value of Q different values of constants T. Th 
ih 
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and also the coefficients of the characteristic polynomial aos ay 


seeeeeee BL have to be calculated, since the initial conditions of 


the synchronous machine vary with different values of Q. 


Example for P = 1.0 and Q= .6 


To calculate Ty Th)» To Ty» Th Ty, the set of equations (30) is used 


and the values are found to be 


le i On 

a 
Tt = Bo L208 5 ue 
Te 16,04 + ous 

® s 


~@O3h. 7a. = D 


¥d mevig hos sotizss beswhorsat sean? ats ated Seay stsh eateve ot 
2 nt ef qmtd baw) coulsy 4uiq al ows exasseoo EEA Ae) 
Pie aiid: Bh hobs bps taieor Sia (69) enobjaups io Soe af galled 
bas 0 20 pails jasvstith Lis 102 Sesetuoleo say Sabdonm evemorie 
gted bosglddua ot 1 x2beeqgh ad mode ie 
‘4? 4,7 einsiendo 20 aides snort? ib 0 tc sti “reyst bb ite 101 
gl mre ons ens gb agit Mahi vhlape lak, Og ee ate “gh ag? | 
pe ene | Isimoavleg ptaatzotoarats: sii to eanstatiases add cate bite 
Yo endisibne> tatiini edd aonte ,bs2atunias ed ot sved 8 Licsesee | 


a 
«0 0 esulev sns1st2lb Adiw viav sntdosm evonurdonye Sad 


8. =O bas Ol =o wt sigue’ 


Beau et (OE) antokteypo to tse 942 QT ,oT 4 Ta T yt o,7 samtvolas at 


SS. = e_7 
38 €8 + &L8 «= a 


41 


Ty =i Gi) 22,5 ut 
Th = 6 
Ty, 25.4.5 


The next step is to calculate constants Ay, Ay eeeceeee Agy Age 


The set of relationships given by (31) will serve this purpose and 


the constants are 


A, © = 4239 A, = .4755 
A, = RBS Ay = = .9244 
A, = .339 A, = - .783 
A, = - .126 A, = 1,058 
A, = - .263 


To calculate the constants B,, B, ..+..-s-- B, the set of equations 


given by (29) is used and the constants are 


Bae 1 6.33 Ba 15.667 
Boer ok Be = 21397 
DS gaat ae) Beogt ae) 


With these constants the coefficients of the characteristic poly- 
nomial for the dynamic system can be calculated. The set of equations 


(28) will serve this purpose, giving coefficients of the characteristic 


polynomial as follows: 
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a tag stl), 

Bye 299 Oe a 

a, = 166.479 + 40.232 OE -0201 UL 
@, (iid 414. debostiel 32 2359 ibs ke Lied us 
a, = 358.61 + 302.427 se ae jes Wig) ue 
a, = 186.874 + 152.042 ce jie aD Us 
eye Si PC cee Eee) fs) We rod, ae 

a, = .263 Ub ot 79 


As was mentioned earlier, it is convenient to leave these coefficients 
in terms of be and ue (the overall stabilizer and regulator gain, 
respectively) in order to draw the stability regions on this parti- 
cular plane. When the coefficients are in this form, the program for 
digital computer shown in Appendix II can be used. This program em- 
ploys the conditions for stability directly in the form given by (24) 
and derived in the text previously. This program checks for stability 
for as many points on wee plane as required and specified 
in the program. All points on Hey plane which satisfy the 
stability conditions constitute the stable region on this plane. 

Using the information given by digital computer, results can be dis- 
played graphically as in Fig. 3a. In this figure the stable region 

of the regulated synchronous machine at P= 1.0 and Q= .6 p.u. 


connected to an infinite bus is shown. All points which lie inside 
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the triangular region satisfy the stability conditions (24) and the 
points outside this region are unstable, i.e. they do not satisfy the 
conditions for asymptotic stability. 

Using this figure for a chosen value of the overall stabilizer gain, 
it is possible to find the maximum allowable voltage regulator gain; 
the minimum allowable voltage rezulator gain shows fairly constant 
value for all possible settings of we 

Using the same procedure, the stability regions for all other power 
factors can be established (Q = .3; 0; -.3). All necessary constants 
and the coefficients of the characteristic polynomial were calculated 
and the stability regions established. The results are graphically 
shown in Fig. 3b and 3c, where the effect of different power factor 

on stability region of dynamic power system can be seen. 

It can be seen in the above mentioned figures that the minimum 
allowable voltage regulator gain is constant in all investigated cases. 
Both the maximum voltage regulator and stabilizer gains decrease their 
values and reach certain minimum and then again increase their values 
as the system goes from lagging power factor (positive Q) towards 


leading power factors. 
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5.3 Effect of Different Damping Coefficients on Stability Region 


In this study the real and reactive power output of the synchronous 


machine is constant during the investigation at 
P = 1,0 Q = .6 


The damping coefficient D is varied and the effect of different 
damping coefficients on stability region studied. During the study 


these three distinct values for damping coefficient are assumed: 
Pa lear. See De We 


All other data introduced by (25) are kept constant during the investi- 
gation. 

The results of this study are graphically shown in Fig. 4, where the 
effect of different damping coefficients on stability region of the 
dynamic system under investigation can be seen. 

From Fig. 4 it can be seen that the larger damping coefficient has a 
favourable effect on voltage regulator and stabilizer gain settings. 
As it is expected, the stability region is larger for larger damping 
coefficient and it allows larger maximum allowable gains. The study 
shows that there is a very distinctive change in both maximum 
allowable (voltage regulator and stabilizer) gain settings if the 


damping coefficient is changed. 
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5.4 Stabilizer Time Constant Effects on the Stability Region 


The real and reactive power output of the machine is constant at the 


following values: 
Ree=aeleQ Q = 0 


The stabilizer time constant is varied during the investigation and 
the effect of this variation on stability region of the dynamic system 
is studied. During this study three distinct values of the stabilizer 


time constant are used: 


T. ol, Deed ie 6) Seconds) 


All other data used in this study are those introduced by (25) and 
they are kept constant during this investigation. 

The results were processed graphically and they are shown in Fig. 5. 
In this figure the stabilizer time constant effect on the stability 
region of the dynamic system under investigation can be seen. 

From Fig. 5 it can be seen that the minimum allowable voltage regula- 
tor gain is constant in all investigated cases. In the high ne 
region, a larger stabilizer time constant allows a smaller voltage 
regulator gain at the low uy portion. The maximum allowable 


stabilizer gain has the same value in all investigated cases. 
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5.5 Exciter Time Constant Effect on the Stability Region 


In this study the real and reactive power output of the machine is 


kept constant at 
P = 1.0 Q = O 


The exciter time constant is varied during tiis investigation and the 
effect of this variation on stability region of the investigated 
dynamic system is studied. The exciter time constant is given the 


following values: 


ie =e 40s Lo (seconds) 


During this investigation all additional data are supplied by (25) 

and they are kept constant. 

The results are shown graphically in Fig. 6. In this figure the 
exciter time constant effect on the stability region of the dynamic 
system under study can be seen. 

In general, the larger the exciter time constant the smaller the 
maximum allowable voltage regulator gain. The minimum voltage regula- 
tor gain shows constant value in all investigated cases and so does 


the maximum allowable stabilizer gain. 
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5.6 One-Time Constant Governor Effect on the Stability Region 


As in the previous study, the real and reactive power output of the 


machine is kept constant at 
Pe) 3170 Q = 0 


The effect of a one-time constant governor is studied here and there- 
fore this constant is varied in this study. The following three 


values are considered: 
Tt =) Of Li 1.4) (seconds) 


Fig. 7 shows the one-time constant governor effect on the stability 
region of the dynamic system under study. 

From Fig. 7 it can be seen that at low ut portion of the high us 
region, a fast-acting governor allows only a smaller voltage regulator 
gain than a slow-acting governor. There is no difference between 
cases in the low region and the maximum allowable stabilizer gain 


shows constant value. 
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5.7 Two-Time Constant Governor Effect on the Stability Region 


The real and reactive power output of the machine is kept constant at 


these values: 
P=) 120 Op tee (8) 


Two-time constant governor effect is studied here and therefore this 
constant is varied during this study. The following three distinct 


values were chosen for the purpose of the study: 


Tee =) O03 L; 2) seconds) 


The results are shown in Fig. 8. This figure shows the two-time 
constant governor effect on the stability region for dynamic system 
under study. From Fig. 8 it can be seen that for a system with a 
two-time constant governor the effect of the values of the time 
constant on the voltage regulator gain settings is smaller than that 
of a system with a one-time constant governor. The minimum voltage 
regulator gain and the maximum allowable stabilizer gain is constant 


and there are only very small changes in the low ut region. 
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CHAPTER VI 


CONCLUSIONS AND SUGGESTIONS 
FOR FURTHER RESEARCH 


Ge Summary and Conclusions 


The dynamic system studied in this thesis consists of a regulated 
synchronous machine connected to an infinite bus. Using the trans- 
formation of the system into Schwarz form (14) and the necessary and 
sufficient conditions of a dynamic system established in this thesis 
(24) the closed stability regions on ur us plane were obtained. 

The method presented here is fairly simple and can be used for both 
designing and power system operational practice. The closed stability 
regions obtained here should be of some interest as no paper has been 
published so far on this subject. 

The method of approach used in this thesis is quite general and gives 
the opportunity to study the effect of many parameters on the stability 
region of a dynamic system. Six parameters were chosen for this 
purpose in this thesis and the studies were conducted and introduced 
earlier in the text. The following general results are observed in 
these studies. 

It can be observed that the general shape of the stability region for 
the system under study is triangular. This is an agreement with the 
operational power system practice, where for practical purposes 
smaller values of the overall stabilizer gain and higher values of 


the overall regulator gain are chosen. There seems to be a definite 
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trend to stay in the lower part of the stability region in order to 
ensure good controllability of a power system. It can be verified 

by [3] or [4] where the voltage regulator data are introduced for the 
systems studied in those papers. 

It can also be seen that in the studies presented in this thesis most 
of the changes in the stability regions are in the low EE portion of 
the high ue region. The minimum allowable voltage regulator gain 
shows fairly constant value in all presented studies. 

Once the stability conditions for the system under study are obtained, 
the rest of the calculations (constants and coefficient of the 
characteristic polynomial) can be computerized. The study as presented 
here does not require special output or computer storage facilities. 
Generally, all results and conclusions presented earlier are in 
agreement with [3]. In that study, however, open stability regions 
were obtained and the system under investigation in that paper had 
different parameters than the dynamic system in this study. Therefore, 


only a general comparison of these two studies is possible. 


6.2 Suggestions for Further Research 


The stability study of the individual unit is presented above 

and necessary and sufficient conditions are established, assuming the 
rest of the system as an infinite bus. These conditions can be com- 
pared with conditions for stability of the system with more than one 
generating unit in order to find any possible variations in both 


stability conditions and stable regions of the system. 
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For practical studies more information on the variation of system 
parameters would be useful and more effective results can be obtained. 
If the limits within which a parameter varies were known,some domain 
of attraction could be established for the investigated system.This 
would be the common area for all stability regions obtained by varying 
the system parameters within the assumed limits.This information can 
be of vital importance in operational practice as well as in design 


of power systems. 
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The real power output of the machine is kept constant in all cases 


(P - 1.0 p.u.) and the reactive power output of the machine (Q) is 


varied, showing the following initial conditions: 
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APPENDIX IL 
FORTRAN IV G PROGRAM 
C PROGRAM TO FIND THE STABILITY REGIONS 


C OVERALL STABILIZER GAIN IS X 
C OVERALL EXCITATION GAIN IS Y 


0001 DO 6 J=1, 32 
0002 DO 7 T=1, 32 
0003 X=(-3+1) 

0004 Y=J 

0005 WRITE(6,10)X,Y 


C COEFFITIENTS OF CHARACTERISTIC POLYNOMIAL WILL FOLLOW 
0006 A=.594 


0007 B=(35.361+(.53*X) ) 
0008 C=(123.648+( 31. 36*X)-(.597*Y) ) 
0009 D=(345.607+(97.667*X)—-(1.218*Y) ) 
0010 E=(306.133+(101. 4*X) -(2.187*Y) ) 
0011 F=(163.277+(135.087*X)-(6.511*Y) ) 
0012 G=(53.15-(1.975*X)-(1.675*Y) ) 
0013 H=(.263*Y)-.79 
C CONDITIONS FOR STABILITY WILL FOLLOW 
0014 CONDI=(B/A) 
0015 WRITE(6,11) CONDI 
0016 COND2=(B*C)-(A*D) 
0017 WRITE (6, 11) COND2 
0018 Bl=((B*C—A*D) / (A*B) ) 
0019 B2=((D*E-C*F) / (A*D) ) 
0020 B3=((F*G—E*H) /(A*F) ) 
0021 Cl=((B1*(D/A) )-(B2*(B/A) )) /B1l 
0022 C2=((B2*(F/A) )-(B3*(D/A) )) /B2 
0023 C3=(H/A) 
0024 Dl=((C1*B2)~(B1*C2))/C1 
0025 D2=((C2*B3)-(B2*C3))/C2 
0026 El=((D1%*C2)-(C1*D2))/D1 
0027 E2=(H/A) 
0028 COND 3=((D/B)-(B2/B1) ) 
0029 WRITE (6, 11) COND3 
0030 COND4=((B2/B1)-(C2/C1)) 
0031 WRITE(6,11) COND4 
0032 COND5=((C2/C1)-(D2/D1) ) 
0033 WRITE(6,11) CONDS 
0034 COND6=((D2/D1)-(E2/E1)) 
0035 WRITE(6, 11) COND6 
0036 COND7=(E2/E1) 
0037 WRITE(6,11)COND7 
0038 10 FORMAT(1H ,2F13.4) 
0039 11 FORMAT(1H ,E14.7) 
0040 7 CONTINUE 
0041 6 CONTINUE 
0042 STOP 


0043 END 
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